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LABORATORY  MEASUREMENTS  OF  SCATTERED  ELECTROMAGNETIC 

RADIATION  FROM  TWO  DIMENSIONAL  METALLIC  AND 
DIELECTRIC  ROUGH  SURFACES. 


ABSTRACT 

Bistatic  and  non  specular  scattering  cross  section  measurements  of  C02  laser  radiation  from  roughened 
metal  and  dielectric  surfaces  were  made  and  compared  to  predictions  given  by  the  facet  (tangent  plane 
approximation)  model.  The  incident  radiation  was  linearly  polarized  perpendicular  to  the  incident 
plane.  The  scattered  polarization  state  was  analyzed  along  two  directions,  perpendicular  (HH)  and 
parallel  (HV)  to  the  scattering  plane.  For  the  first  time,  such  polarization  dependent  bistatic 
measurements  are  reported  over  the  entire  hemisphere,  using  a  scattering  apparatus  designed  especially 
for  this  work.  The  average  slopes  and  radii  of  curvature  of  the  roughened  surfaces  utilized  in  this 
study  were  determined  from  surface  tracings  using  a  mechanical  profilometer.  The  (Oq^HH  nulls 
predicted  by  the  facet  model  have  been  verified  at  steep  incident  angles  for  both  metallic  and  dielectric 
surfaces.  Agreement  is  generally  better  for  metals  than  dielectrics.  At  shallower  angles  measurements 
diverge  from  theoretical  predictions.  Our  data  indicate  that  the  departure  from  the  predictions  of  the 
facet  model  is  most  likely  associated  with  surface  curvature.  A  number  of  calculations  of  polarization 
dependent  scattering  cross  sections  from  metal  and  dielectric  cylinders  of  radius  r  as  a  function  of  r (k 
have  been  performed.  These  calculations  suggest  that  even  for  large  radii  surfaces,  rapid  amplitude 
and  phase  variations  are  responsible  for  the  onset  of  depolarization  at  shallower  incident  and  scattering 
angles,  leading  to  the  disagreement  with  the  tangent  plane  model  predictions.  They  also  clearly 
demonstrate  why  the  facet  model  is  a  better  approximation,  in  the  region  of  validity,  for  metals  than 
dielectrics.  The  role  of  curvature  in  determining  the  operable  regions  of  the  tangent  plane  model  is 
further  illuminated  through  a  series  of  scattering  measurements  from  metal  wires  with  radius  of 
curvature  on  the  order  of  X.  The  experimental  measurements  are  described  in  detail  and  data  for  both 
roughened  metals  and  dielectrics  are  presented  for  several  scattering  aspect  and  depression  angles. 


1 


INTRODUCTION 


✓ 

Scattering  of  electromagnetic  radiation  by  rough  surfaces  has  been  the  subject  of  great  interest 
for  several  decades!1!.  Interest  in  this  derives  from  a  need  to  study  terrain  characteristics  of 
inaccessible  sites  and  to  detect  and  characterize  small  deviations  from  optically  smooth  surfaces.  In 
either  case,  the  goals  are  to  relate  electromagnetic  scattering  data  to  the  properties  of  the  scattering 
surface.  In  principle,  knowledge  of  the  exact  shape  of  the  scattering  surface  and  the  total  electric  field 
E  at  the  surface  allows  one  to  compute  the  electric  field  E2  at  the  point  of  observation.  The  exact 
solution  is  given  by  equation  (1)  where  k  is  the  magnitude  of  the  wave  vector  of  the  incident  radiation 
and  R  is  the  distance  between  some  point  on  the  illuminated  surface  and  the  observation  point!2! : 


where 


E2(P)=i{(E^-'^>dS 


(1) 

(2) 


A  cursory  look  at  equation  (1)  immediately  confronts  one  with  the  two  fundamental  obstacles 
to  finding  E2(P).  One  problem  is  that  the  exact  shape  of  the  surface  is  unknown.  Secondly,  given  the 
detailed  geometry  of  the  surface,  there  still  remains  the  problem  of  knowing  the  value  of  the  total  E  on 
the  surface,  as  given  by  (3). 


E  =  E|  +  E,  (3) 

Ei  and  E,  represent  the  incident  and  scattered  electric  fields  on  the  surface ,  respectively.  To 
obtain  the  value  of  E,  one  must  solve  the  boundary  conditions  for  both  E  and  H  at  the  surface. 
Unfortunately,  closed  form  solutions  for  the  latter  problem  exist  only  for  plane  surfaces.  In  the  high 
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frequency  or  geometrical  optics  (GO)  limit,  the  ratio  of  scattered  to  incident  electric  field  is  the  same  as 

% 

the  corresponding  quantities  obtained  from  the  Fresnel  coefficients.  Consequently,  in  the  GO  limit, 
this  part  of  the  problem  is  eliminated.  The  need  to  characterize  the  surface,  however,  remains.  An 
exact  mapping  of  the  surface  is  an  insurmountable  task  and  may  not  be  necessary.  It  is  usually 
assumed  that  different  microscopic  surface  shapes  will  yield  similar  scattering  patterns,  when 
illuminated  over  a  sufficiently  large  area  of  the  rough  surface.  The  statistical  characterization  of  rough 
surfaces  is  motivated  by  the  need  to  obtain  closed  form  expressions  for  E2(P)  without  detailed 
knowledge  of  the  surface  shape.  A  model  in  which  one  assumes  a  random  distribution  of 
hemispherical  bosses!3!,  the  tangent  plane  model  (TP)f4l  and  the  two  scale  roughness  model!5!,  are 
commonly  used.  The  TP  approximation  is  the  most  straightforward  approach.  In  the  GO  limit,  all 
shapes  can  be  handled  by  the  TP  approximation  in  which  only  the  average  surface  slopes  are  of 
interest.  To  obtain  a  statistical  representation  of  the  surface,  it  is  important  that  the  illuminated  area 
provides  an  adequate  representation  of  the  surface  characteristics.  Papa,  Lennon,  and  Taylori6!, 
henceforth  denoted  as  PLT,  have  shown  that  for  uncorrelated  heights  and  slope  distributions,  in  a 
Gaussian  surface,  the  physical  optics  approach  yields  identical  results  to  those  obtained  in  the  GO 
limit,  provided  that  the  correlation  length  T  is  much  larger  than  the  wavelength  and  the  average  slopes 
are  significantly  less  than  one.  Another  important  observation  of  PLT  is  that  in  the  TP  regime,  the 
polarization  of  the  scattered  wave  is  independent  of  the  detailed  statistical  properties  of  the  surface. 


Previous  measurements  of  rough  surface  scattering!7*8!  were  performed  in  the  monostatic 
configuration  or  under  conditions  where  the  scattered  radiation  was  in  the  plane  of  incidence.  Such  an 
arrangement  precludes  the  experimental  study  of  some  interesting  features  of  polarization  dependence. 
These  manifest  themselves  only  in  configurations  where  the  scattering  plane  is  different  from  the 
incident  plane.  The  apparatus  utilized  in  this  study  allows  the  measurement  of  scattered  radiation  in 
both  the  polar  and  azimuthal  directions. 
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Unlike  smooth  surfaces,  roughened  surfaces  scatter  in  all  directions.  This  requires  the  detector 
to  be  moveable  over  a  hemisphere.  Practical  considerations,  especially  for  a  liquid  nitrogen  cooled  IR 
(infrared)  detector,  limit  the  motion  of  the  detector  to  the  horizontal  plane  in  the  laboratory.  To  insure 
accessibility  of  arbitrary  scattering  angles,  the  polar  angle  8  s  and  particularly  <^,  the  azimuthal  angle, 
require  that  the  target  surface  be  free  to  rotate  around  two  perpendicular  axes.  At  the  same  time,  the 
polarization  stare  of  the  transmitted  radiation  must  be  adjusted  to  maintain  a  well  defined  polarization 
with  respect  to  the  surface  normal.  A  description  of  the  scattering  system  and  the  correspondence 
between  laboratory  angles  in  which  the  mean  surface  normal  changes  direction,  and  the  conventionally 
defined  field  angles,  is  presented  m  section  IL 

An  algorithm  to  compute  laboratory  angles  for  a  specific  set  of  field  angles  is  presented  in  the 
Appendix.  Data  were  collected  in  a  bistatic  configuration  where  the  incident  and  scattered  planes  were 
different.  The  co-pol  and  cross-pol  scattering  cross  section  measurements  were  compared  with  the 
theoretical  predictions  of  PLT^6L  Using  a  form  of  the  scattering  cross  section  given  by  Barrick!4! 
(equation  4),  PLT  describe  the  angular  dependences  of  o0  as  a  function  of  polarization,  surface 
roughness,  and  dielectric  constant  The  scattering  cross  section  is  given  by!4-6!; 

So  =|  Ppq|  JS  (4) 


In  equation  (4),  S  is  the  shadowing  function!4*6!  and  J  is  the  probability  density  function!4-6!  for 

_  ^ 

surface  slopes.  J  is  proportional  to  (-*-)  times  an  exponential  function  of  ,  where  T  is  the  average 

O’  o  ° 

facet  spacing  and  a  is  the  average  facet  depth. 

As  described  by  Bamck  ,  a  specular  point  in  the  surface  reflects  like  a  tilted  plane  tangent  to 
the  surface  at  that  point  The  (Jpq  in  equation  (4)  are  given  below  in  terms  of  the  Fresnel  coefficients 
for  HH  and  W  waves  and  spherical  scattering  coordinates  (see  Figure  2). 
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q  a2a3Rn(i)  +  sm(0i)sin(0I)sin2(0I)R1 

Pw  *  aia4 

(5) 

a  sin  0,[-sin  0ia3R»(i)  +  sin  0ta2Rj.(i)] 

PHV  *  aia4 

(6) 

_  sin  0,[sin  0sa2Rii(i)  +  sin  O^R ,  (i)  ] 

PvnH  =  aja4 

(7) 

_  -sin  0jsin  0,sin2  4»«Rii(i)  -  a2a3R  •  (i) 

Prat  -  aia4 

(8) 

with 


ai=l  +  sin  0jsin  0,cos  0,  -cos  0iCos  0f 

(9) 

a2=cos  0jsin  0t  +  sin  9;cos  0fcos  <J>, 

(10) 

a3=sin  Ojcos  0,  +  cos  0jsin  0,cos  0, 

(ID 

a4=cos  0-,  +  cos  0, 

(12) 

The  angle  (i)  is  the  angle  of  incidence  with  respect  to  the  local  normal  of  the  facet  It  is  defined 
below  in  terms  of  the  scattering  angles  of  Figure  2. 


cos  i  =  ^{1  -  sin  0;sin  0,cos  0,  +  cos  0jcos  0S]3 


(13) 


In  equation  (4),  the  J  terni  is  proportional  to  the  average  number  of  facets  having  slopes  which 
scatter  into  the  observation  direction.  The  shadowing  function  S  gives  the  fraction  of  the  total  number 
of  specular  points  not  shadowed!9!.  In  the  analysis  of  PLT,  it  is  pointed  out  that  neither  the  shadowing 
function  nor  the  slope  statistical  J  term  has  an  influence  on  the  position  of  nulls  in  the  scattering  cross 
section.  Thus,  the  predicted  angular  position  of  nulls  in  a0  as  given  by  these  authors  is  based  solely 

on  the  behavior  of  |  Hpqp  in  equation  (4).  Furthermore,  their  analysis  shows  that  for  a  given  incident 
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signal  polarization,  there  exists  nulls  in  the  co-pol  scattering  cross  section  for  various  angles. 

/ 

.  * 

In  this  work,  the  predicted  polarization  dependent  behavior  of  the  scattering  cross  section  as 
given  above  was  observed,  for  a  subset  of  incident  scattering  angles.  For  steep  incident  and  receive 
polar  angles,  the  data  given  in  this  report  conforms  to  predictions  of  the  TP  approximation. 
Deficiencies  in  the  TP  model  appear  for  larger  angles.  The  inability  of  the  TP  approximation  to 
accurately  describe  the  scattering  of  radiation  for  all  angles  i  has  been  noted  previously  in  other, 
monostatic  measurements!7*8!.  The  departure  from  the  TP  model  has  been  attributed  to  the  neglect  of 
small  scale  roughness  of  the  surface  and  to  neglect  of  multiple  scattering  from  facets  with  large  slopes. 
A  convincing  argument  can  be  made  that  curvature  inherent  in  rough  surfaces  is  a  major  cause  of  the 
observed  deviation  from  the  facet  model. 

A  good  way  to  test  the  adequacy  of  the  TP  approximation  is  to  juxtapose  the  TP  predictions 
with  exact  calculations.  For  example,  we  can  compare  the  reflection  coefficients  computed  for  the  case 
of  radiation  incident  on  a  cylinder  with  the  results  obtained  if  we  approximate  the  cylinder  by  a  set  of 
tangent  planes.  It  is  possible  to  describe  analytically  the  polarization  dependent  scattering  amplitude 
and  phase  for  perfect  curved  geometries,  such  as  cylinders,  in  terms  of  the  parameter,  r/X.  The  onset 
of  depolarization  for  certain  scattering  geometries  and  the  departure  of  data  from  the  TP  model  is 
shown  to  be  critically  dependent  on  this  parameter.  A  scries  of  calculations  of  scattering  from  metal 
and  dielectric  cylinders  is  given  in  section  IV.  These  calculations  disclose  the  very  rapid  phase  and 
amplitude  variation  of  the  W  and  HH  components  of  the  scattered  field  as  a  function  of  incident  angle 
and  surface  curvature.  A  description  of  these  calculations  and  their  significance  in  understanding 
scattering  from  rough  surfaces  is  also  contained  in  section  IV. 

A  degree  of  insight  into  polarization  dependent  scattering  can  be  attained  from  these 
calculations  by  characterizing  the  roughened  surfaces  according  to  slope  and  radius  of  curvature. 
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Metal  and  dielectric  rough  surfaces  utilized  in  this  study  were  prepared  in  the  following 
manner.  The  surfaces  were  blasted  with  glass  beads  under  a  pressure  of  approximately  80  psi.  The 
bead  diameters  ranged  from  50  to  500  jim.  Profilometer  measurements  were  performed  on  both  metal 
and  dielectric  rough  surfaces.  A  digitization  technique  in  which  the  plotted  profilometer  data  was 
converted  into  computer  readable  data  files  was  employed  and  programming  was  developed  to 
calculate  the  average  surface  slope  and  radius  of  curvature.  Utilizing  these  techniques,  the  average 
values  of  these  surface  parameters  for  the  entire  surface  were  determined  as  well  as  that  percentage  of 
the  surface  falling  within  a  certain  range  of  the  parameters.  Table  I  displays  results  of  such 
measurements  for  the  surface  utilized  in  th;s  study  in  which  the  radii  of  curvature  are  given  in  units  of 
the  wavelength.  Finally,  results  from  randomly  arranged  metal  wires  with  r  ~  X  are  reported,  yielding 
further  evidence  that  the  departure  from  the  facet  model  predictions  is  due  to  surface  curvature. 

The  organization  of  this  report  is  as  follows. 

Section  II  contains  a  discussion  of  the  laboratory  setup.  Section  III  is  devoted  to  a  set  of 
laboratory  scattering  data.  In  Section  IV  we  discuss  the  data. 

Appendix  A  contains  formulae  for  transforming  field  angles  into  laboratory  angles.  Appendix 
B  contains  additional  scattering  data. 
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IL  EXPERIMENTAL  SETUP  AND  DISCUSSION  OF  LABORATORY  ANGLES 


LN2 

HgCdT*  MnK 


rough  •urfaet 

Figure  1  Component  Diagram  of  Scatterometer  System 
The  source  of  radiation  is  a  C02  laser  (Ultra  Lasertech  Model  5122)  providing  10  watts  of 

linearly  polarized  radiation  with  a  choice  of  48  lines  between  9.2  to  10.8  microns.  The  ability  to  rotate 
the  linearly  polarized  radiation  to  any  desired  angle  is  provided  by  a  K-rotator.  The  rough  surface  is 
allowed  to  rotate  through  two  degrees  of  freedom  which,  in  conjunction  with  the  position  of  the 
detector,  simulates  the  environment  of  a  bistatic  radar  in  the  field.  At  the  receive  end,  discrimination 
between  polarization  states  is  accomplished  by  rotating  a  Brewster  plate  analyzer  and  wire  grid 
polarizer  analyzer.  The  combination  of  the  two  yields  an  extinction  ratio  of  better  than  1000  :  1.  The 
radiation  is  collected  by  a  two  inch  diameter  ZnSe  lens  with  a  focal  length  of  20  inches.  A  liquid 
nitrogen  cooled  HgCdTe  detector  (detectivity  =  lO1^  cm  Hz^/Watt),  mounted  on  a  rotatable  detector 
arm  in  the  focal  plane  of  the  lens,  scans  the  resultant  intensity  pattern.  Scanning  over  a  finite  angular 
spread  and  averaging  the  data  is  necessitated  by  the  scintillation  pattern  that  results  from  the  narrow¬ 
band  laser  radiation  scattering  from  the  rough  surface.  The  detector  output  is  fed  to  a  lock-in  amplifier 
and  the  angle  averaged  signal  displayed  on  a  Mac  II  computer. 
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The  rough  surface  is  mounted  on  a  goniometer  which  can  rotate  through  a  range  of  90° .  This 
rotation  is  about  an  axis  that  is  parallel  to  a  horizontal  surface.  The  goniometer  is  mounted  on  a  turn¬ 
table  which  rotates  through  a  vertical  axis.  By  means  of  these  two  components,  the  target  surface  can 
be  rotated  around  two  axes. 

The  detector  has  a  range  of  270°  and  is  rotated  independently  of  the  target  surface. 

The  setup  can  simulate  the  three  independent  field  angles,  ie,  the  incident  and  scattering  polar 
angles  from  0  to  90  °,  and  the  scattering  azimuthal  angle  from  0  to  180°.  The  rest  of  the  azimuthal 
range,  180-360°,  can  be  obtained  by  rotadng  the  target  surface  around  its  mean  normal. 

The  transformation  of  the  field  angles  into  the  corresponding  laboratory  angles  will  now  be 
presented. 

In  the  field  frame,  the  target  surface  orientation  is  fixed  and  the  transmit  and  receive  directions 
(vectors)  can  vary  independently  over  a  hemisphere  (see  Figure  2).  In  this  frame,  we  designate  the 
following  angles  0j,  0S,  and  <j>s.  0j  stands  for  the  incident  angle,  the  angle  between  the  incident 
propagation  direction  and  target  normal. 


H  •  n  =  -Cos(0O  (14) 

where  r.  is  a  unit  vector  along  the  direction  of  propagation,  and  n  is  a  unit  vector  along  the  surface 
normal.  Similarly, 


r,  •  n  *  Cos(0,) 


(15) 


9 


where  r  is  a  unit  vector  along  the  detector  direction  and 
s 


[£j '  "<ri  •  ■)!  .  tfj  -  n(r»j  jOl  = 

Sin(6i)  Sin(6t) 


(16) 


10 


In  the  laboratory  frame  (Figure  3),  the  propagation  direction  is  fixed,  which  is  taken  along 
the  negative  z-axis.  Thus, 


fi(L)  =  -i3 


(17) 


where  L  designates  the  laboratory  frame  and  i3  is  a  unit  vector  along  the  z-axis  (Figure  3).  The 
detector  arm  is  in  the  x-z  plane  and  can  be  rotated  through  270°  (Figure  4). 


Figure  4.  Laboratory  system  for  \j/  =  0  and  2;  =  0  showing  the  detector  position  rD(L) 

The  target  surface  normal  can  be  rotated  around  two  perpendicular  axes.  One  rotation  is 
around  the  y-axis  (Figure  4),  the  axis  perpendicular  to  the  x-z  plane  containing  the  transmitter  and 
receiver.  The  angle  describing  this  rotation  is  designated  by  !;.  The  range  of  2;  is  0°  to  90°.  The 
detector  arm  requires  a  180°  range  for  any  fixed  2;.  Hence,  the  full  range  of  the  detector  has  to  be  270°. 
Figure  5  shows  the  case  for  ^  *  0.  The  relation  between  2;s  and  T|  is  £s  =  “n  -  90°.  2;  equals  the 
incident  field  depression  angle,  0},  when  the  target  is  rotated  around  the  y-axis  only. 
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A  second  rotation  can  be  performed  around  the  line  which  is  the  intersection  of  the  target  plane 
and  the  x-z  plane  (Figure  3). 

The  first  rotation  fixes  \  and  the  second  rotation  fixes  y.  When  the  target  surface  is  rotated 
through  ^  and  \p  in  succession,  the  normal  to  the  target  surface,  expressed  in  the  unit  vectors  of  the 
laboratory  coordinate  system  is: 

n"  =  -ii  Sin(^)  Cos(v)  +  h  Sin(y)  +  i3  Cos  (£)  Cos(y)  (19) 

While  4i  is  now  not  equal  to  the  incident  depression  angle,  we  will,  for  convenience,  refer  to  it  as  the 
"laboratory  depression  angle",  and  similarly  refer  to  y  as  the  laboratory  aspect  angle. 


*2 

A 


Figure  5.  Laboratory  system  for  ^  *  0  and  y  =  0 
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We  note  the  scalar  products  (14)-(16)  defining  die  various  field  angles  retain  die  same  form 
when  n*  is  substituted  for  n.  The  following  expressions  relate  laboratory  angles  to  field 

angles  (0j,  0S  ,4>s) 

Cos(0i)  *  -r*,  •  n"  =  Cos  (£)  Cos(y)  (20) 

Cos(0,)  =  rDn"  =  Sin(rj-^)  Cos(y)  (21) 


and 


„  v  fi-fD  -  (rs-  n'Xro-n') 
Cos(4>5)  =  ==  7  ==  • 

V 1  -  (ri*  n*)2  V 1  -  (ro-  iO2 
Substituting  Eqns.  (17),(18),(20)  and  (21)  into  Eqn.  (22)  yields: 


(22) 


Cos(<t>5)  = 


-SinCn)  +  Cos2(y)  Sin(^-n)  Cos(^) 

V  l-Cos2(£)  Cos2(\|f)  V  l-Cos2(v)  Sin2(^-T|) 


(23) 


An  algorithm  for  computing  laboratory  angles,  T|,  \p  for  a  given  set  of  field  angles,  0j,  0S> 
4>S'  can  be  found  in  the  Appendix. 
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m  Scattering  Data 


Scattering  data  were  obtained  from  roughened  aluminum  surfaces,  roughened  dielectric 
surfaces,  randomly  arranged  wires,  and  a  randomly  distributed  set  of  glass  beads. 

A  sand  blaster  was  used  to  roughen  aluminum  surfaces.  The  abrasive  material  used  in 
preparing  the  surfaces  was  a  collection  of  glass  beads  ranging  from  SO  to  500  microns  in  diameter. 
The  results  of  profilometer  tracings  indicate  that  the  roughening  produces  randomly  distributed  valleys 
ranging  from  3  to  9  microns  with  an  average  spacing  of  70  to  100  microns.  Examination  of  the 
surfaces  under  a  microscope  indicated  the  existence  of  sharp  edged  patches.  To  remove  these  sharp 
edges,  the  surfaces  were  electropolished,  removing  approximately  1  micron  of  material  from  the  upper 
surface.  Statistical  data  is  given  in  Table  I  indicating  radius  of  curvature,  slope,  and  average  facet 
depth  for  the  surface  (EBT  4+13)  studied.  The  histogram  characterizing  this  surface  according  to 
distribution  of  facet  slopes  is  plotted  in  Figure  6. 

Data  are  presented  for  both  copol,  o0(HH),  and  cross  pol,  o0(HV),  measurements,  as  well  as 
ratios  of  o0(HH)/a0(HV).  The  results  are  displayed  in  both  laboratory  and  field  coordinates.  In  the 
lab  coordinates,  the  transmit  and  receive  laboratory  depression  angles,  and  £,s  are  set  equal  and 
remain  fixed.  The  target  surface  is  rotated  around  a  horizontal  axis.  This  rotation  defines  the 
laboratory  aspect  angle.  Changing  the  laboratory  aspect  angle  induces  changes  in  both  0j  and  8S>  the 
incident  and  scattered  polar  angles.  However,  there  is  an  intrinsic  advantage  to  performing 
measurements  in  the  laboratory  coordinate  system  since  the  laboratory  aspect  angle  is  identical  to  the 
angle  subtended  by  the  global  normal  and  the  facet  normal  ( recall  that  according  to  the  TP  model,  the 
scattering  at  a  given  angle  is  due  to  facets  for  which  the  specular  condition  is  fulfilled).  Figures  7-8 
display  both  ratios  of  o0(HH)/o0(HV)  and  o0(HH)  and  c0(HV)  in  laboratory  coordinates.  The  data 
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points  plotted  for  o0(HH)  and  o0(HV)  have  been  radoed  in  such  a  way  that  the  data  at  some  small 
aspect  angle  coincides  with  the  model  prediction.  The  average  slope  parameter,  T/o^,  which  enters 
into  the  TP  model  is  obtained  from  the  results  of  profilometer  measurements  (Table  I).  Figures  9-10 
show  data  for  a0(HH)/<J0(HV)  and  a0(HH)  and  c0(HV)  for  configurations  in  which  the  field  angles  0j 
and  0S  and  the  azimuthal  scattering  angle  <J>S  is  varied. 

Figures  11-12  display  scattering  data  for  roughened  plastic  surfaces  with  n  =  1.6  and 

k  *  0.002. 

Roughening  of  the  plastic  surfaces  was  accomplished  in  the  following  manner.  The  plastic 
surfaces  were  softened  in  acetone  and  compressed  between  roughened  aluminum  surfaces  using  c- 
clamps.  After  hardening,  the  roughened  plastic  surfaces  were  measured  using  a  mechanical 
profilometer  and  found  to  have  surface  depths  and  slopes  very  similar  to  the  metal  surfaces  used  in 
their  preparation. 

Figure  13  displays  scattering  data  from  randomly  distributed  glass  beads.  Figure  14  displays 
scattering  data  from  a  glass  surface  roughened  with  glass  beads.  Figures  15-16  show  scattering  data 
from  aluminum  wires  randomly  distributed  on  an  absorbing  flat  surface. 

To  obtain  a  better  understanding  of  both  the  success  and  failure  of  the  tangent  plane  approximation,  we  prepared  a 
sample  target  surface  consisting  of  randomly  distributed  aluminum  wires.  Wire  from  a  spool  was  continuously  wound  two  or 
three  layers  deep  in  a  random  fashion  on  aflat  absorbing  plate.  The  radius  -a-  of  the  wire  was  125  pm  and  the  corresponding 

A/a  was  0.85.  For  such  a  large  value  of  A/a,  the  tangent  plane  approximation  should  not  work  at  all. 
Interestingly,  however,  the  data  for  the  ratio  a0(HH)/c0(HV)  at  a  60°  laboratory  depression  angle  is 
remarkably  similar  to  that  obtained  for  the  randomly  roughened  surfaces.  A  plausible  reason  for  this 
similarity  will  be  presented  in  the  last  section  of  this  report 
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Table  I 


Surface  Slopes  Depths  and  Radii 
of  Curvature  for  EBT  4+13 


Glass  Bead  (Dia.)  j 

■M 50  -  500  microns 

1 

Average  Slope 

8B|| 

0.081 

RMS  Slope 

Wmm 

0.118 

Standard  Deviation 

A£| 

0.086 

— 

Average  Facet  Depth 

BBI 

188  microns 

.  - 

IKBSEB 

HsSl _ 

<0.05 

518 

<0.10 

ftllllg 

793 

<0.20 

95.4 

<030 

98.7 

<0.40 

vx.sXvxNttwx? 

|  993 

<030 

99.9 

%  of  Surface  with  Radius 

of  Curvature  in  Wavelengths 

(microns) 

•  <  | 

<3 

1  3.6 

<4 

1  53 

<5 

IIIIH 

1  7.8 

<6 

103 
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Figure  6  Distribution  of  Surface  Slopes  for  EBT4+13 


Aluminum  data  (Plots  are  in  Laboratory  Angles) 


EBT4+U 


Figure  7a  cr0(HH)/c0(HV)  for  laboratory  depression  angles  -  £s  = 


Figure  7b  o0(HH)/o0(HY)  for  laboratory  depression  angles  ^  = 


HH/HV  ratio  C  HH/HV  rail* 


EBT4+13 


re  7c  ct0(HH)/o0(HV)  for  laboratory  depression  angles  =  £  =40 

1  s 


Figure  7d  a0(HH)/c0(HV)  for  laboratory  depression  angles  £.  =  =  60 
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Ilinilll 


Figure  8a  cr0(HH)  and  o0(HV)  for  laboratory  depression  angles  L  ■  £  *  10 

I  S 


20/20  HH  mod  HV  for  EBT4+13 


■oooi  iiii  m*rTT-n  »  '  »nr*t  i  i  **— r-i  i  i  i  i  i  i  m 
0  $  10  19  20  29  30 

Aspect  angle  (degree*) 


Figure  8b  o0(HH)  and  a„(HV)  for  laboratory  depression  angles  =  Ss  =  20 
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Figure  8c  o0(HH)  and  a0(HV)  for  laboratory  depression  angles  *  E  «  40 

1 


Aspect  mi  git  (degrees) 

Figure  8d  ct0(HH)  and  a0(HV)  for  laboratory  depression  angles  =  4S  =  >0 
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Aluminum  data  (Plots  are  in  Field  Angles) 


Figure  9a  a0(HH)/o0(HV)  for  field  polar  angles  8j  =  9S  =  10 


Figure  9b  co(HH)/c0(HV)  for  field  polar  angles  9;  =  0S  =  20 


Figure  9c  ao(HH)/<r0(HV)  for  field  polar  angles  0,  =  0S  =  40 


Figure  9d  a0(HH)/c0(HV)  for  field  polar  angles  0,  =  0S  =  60 


AZIMUTHAL  ANGLE  (DEGREES) 


Figure  10a  o0(HH)  and  o0(HV)  for  field  polar  angles  Q\  *  0S  *  10 


AZIMUTHAL  ANGLE  (DEGREES) 


Figure  10b  a0(HH)  and  c0(HV)  for  field  polar  angles  8,  =  8S  =  20 
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»  40  DEGREE  POLAR  ANGLE 
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Figure  10c  oc(HH)  and  o0(HV)  for  field  polar  angles  0j  =  0S  =  40 
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Figure  lOd  o0(HH)  and  a0(HV)  for  field  polar  angles  0j  =  0S  =  60 
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Aximtbal  ingk  (dcgracs) 


Figure  11  a0(HH)/c0(HV)  for  field  polar  angles  *  6S  =  7.5 


Figure  12  a0(HH)  and  o0(HV)  for  field  polar  angles  0,  =  0S  =  7.5 


Glass  data  (Plots  are  in  Laboratory  Angles) 


✓ 


Figure  13a  a0(HH)/a0(HV)  for  laboratory  depression  angles  =  £s  =  20 


Figure  13b  o0(HH)/o0(HV)  for  laboratory  depression  angles  ^  =  40 
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hh/hv 
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Ratio  HH/HV 

Glass  5  mlnutaa  abrasion  with  BT4 
doprossion  ■  20/20 


Figure  14a  o0(HH)/a0(HV)  for  laboratory  depression  angles  ^  ^  as  20 


Figure  14b  o0(HH)/o0(HV)  for  laboratory  depression  angles  ^.  =  =  60 


Cron  section  ratio  (HH/ffV)  £=  Cron  taction  ratio  (HH/HV) 


Aluminum  Wire  data  (Plots  are  in  Laboratory  Angles) 


e  15  a0(HH)/a0(HV)  for  laboratory  depression  angles  =  0S  =  20 


Figure  16  o0(HH)/c0(HV)  for  laboratory  depression  angles  0j  =  0S  =  60 
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IV  Discussion  of  Results 


For  the  sake  of  clarity  we  divide  our  discussion  into  four  parts.  We  start  with  scattering  data 
from  metals  and  divide  these  into  two  groups;  data  in  laboratory  angles.  Figures  7*8,  and  data  in  field 
angles.  Figures  9-10. 

As  previously  noted,  the  facet  model  implies  that  the  laboratory  aspect  angle  is  identical  with 
the  angle  subtended  by  the  facet  normal  with  the  surface  (global)  normal.  The  facets  that  radiate  into 
the  detector  all  have  normals  in  the  horizontal  plane,  the  plane  defined  by  the  transmit  and  receive 
directions.  Figures  7-8  show  that  agreement  with  the  TP  model  is  excellent  up  to  15°  aspect  angle  and 
the  onset  of  significant  deviation  from  the  TP  model  occurs  at  an  aspect  angle  somewhere  between  15° 
and  20°.  The  exception  to  this  statement  is  the  deviation  of  o0(HV)  and  the  ratio  c0(HH)/a0(HV)  at 
very  small  aspect  angles.  According  to  the  TP  model,  cr0(HV)  -*  0  as  the  aspect  angle  -4  0. 
Tentatively,  this  deviation  can  be  attributed  to  a  small  amount  (a  few  percent)  of  diffuse  and 
depolarized  scattering  leading  to  a  non  zero  C0(HV)  at  all  angles  including  small  aspect  angles.  A 
plausible  explanation  for  the  origin  of  this  behavior  will  be  presented  later. 

Figures  9-10  exhibit  data  in  field  coordinates  for  fixed  and  identical  transmit  and  receive  polar 
angles  as  a  function  of  the  azimuthal  angle.  Zero  azimuthal  angle  thus  corresponds  to  the  specular 
configuration.  In  the  field  representation,  the  deviation  from  the  TP  model  predictions  occurs  at 
different  azimuthal  angles.  Of  course,  the  underlying  source  of  the  discrepancy  between  the  TP  model 
and  data  must  be  the  same  in  both  laboratory  and  field  coordinate  systems.  By  plotting  the  angle 
between  the  global  and  local  facet  normal  as  a  function  of  azimuthal  angle,  a  set  of  curves  is  generated 
in  which  each  curve  corresponds  to  a  given  polar  angle.  Looking  at  the  data,  one  can  find  the 
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azimuthal  angle  at  which  there  is  a  significant  departure  between  the  TP  model  and  experimental  data. 

These  points  are  plotted  as  x’s  in  Figure  17.  The  x's  lie  approximately  on  a  straight  line  and  show  that 

/ 

in  field  coordinates,  the  departure  of  the  data  from  the  predictions  of  the  TP  model  occurs  when  the 
facet  normal  subtends  an  angle  of  approximately  20°  with  the  global  normal. 


AZIMUTHAL  ANGLE  (DEGREES) 

Figure  17  Facet  Slope  (angle  of  local  facet  normal)  where  Field  Data 
Departs  from  Tangent  Plane  Model 


It  has  been  noted  before  that  there  are  limits  to  the  validity  of  the  TP  approximation  due  to 
curvature.  BrekhovskikhUO]  derived  the  following  criterion  for  the  validity  of  the  TP  approximation: 

47Q-CCOS  $»\  (24) 

where  rc  is  the  radius  of  curvature  and  d  is  the  local  angle  of  incidence.  Unfortunately,  the  inequality 
is  not  specific  enough  to  be  of  use  in  data  analysis. 


We  present  here  a  somewhat  different  argument  In  a  rough  surface  where  the  average  slope  is 
small  (on  the  order  of  0.1  as  is  the  case  here),  larger  than  average  slopes  can  occur  in  several  ways:  the 
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depth  is  larger  than  average  and  the  spacing  between  neighboring  high  and  low  points  has  an  average 
value,  the  depth  is  average  and  the  spacing  between  high  and  low  points  is  smaller  than  average,  or 
finally,  both  depths  and  spacings  are  significantly  different  from  the  average  value.  We  consider  the 
case  indicated  in  Figure  18  below  where  the  depth  is  average  and  the  spacing  between  high  and  low 
points  is  smaller  than  average. 


PI 


Let  the  circle  be  centered  at  a  point  inside  the  surface  that  is  at  2>X  beneath  a  high  point  on  the 
surface  such  that: 


d  =  V  R2  -  (  R  -  h  )2  *  V  2Rh  -  h2  (25) 

If  we  set  R  =  3X  **  31.8  pm  and  h  =  5.76  pm  we  obtain: 

d  =  18.24  pm 

and  0  -  17°. 

Thus,  approximately  half  of  the  facets  corresponding  to  slope  angles  larger  than  17°  will  have  a  radius 
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of  curvature  smaller  than  3X,  the  onset  of  the  resonance  regime!1  ^  for  curved  surface  (cylinder, 
sphere)  scattering  according  to  the  Mie  theory!1  xl 

Surface  characteristics  such  as  depth,  facet  spacing  and  radius  of  curvature,  were  obtained  by 
means  of  a  mechanical  profilometer  "Mitutoyo  Surftest  201",  with  a  5  pm  diameter  diamond  stylus. 
Some  of  the  parameters,  such  as  the  diameter  of  the  glass  bead  and  the  duration  and  pressure  of  the 
blast,  were  varied  during  the  sand  blasting  operation.  Four  differently  prepared  surfaces  were  used  in 
this  study!  13K  After  sand  blasting,  the  surfaces  were  examined  under  a  microscope  and  observed  to 
display  many  sharp  points  and  edges.  Consequently,  the  samples  were  electropolished  to  remove 
these  sharp  features  as  has  been  the  practice  of  other  researchers  in  the  field!12!.  The  surface  data 
given  in  Table  I  is  tabulated  according  to  bead  size  (BT  #).  The  data  displayed  in  Figures  7-8  were 
obtained  for  sample  "EBT  4+13",  an  electropolished  sample  blasted  with  a  mix  of  BT4  and  BT13 
beads.  According  to  Table  I,  the  average  facet  depth  for  this  surface  was  2.88  pm.  This  average  was 
measured  from  the  mean  surface  level  The  quantity  h  that  enters  into  Eqn.(25)  is  double  the  value  of 
the  mean  depth  found  in  Table  I.  Table  I  also  gives  information  about  the  radius  of  curvature  of  facets 
for  all  surfaces  studied.  Radius  of  curvature  calculations  were  performed  by  fitting  a  parabola  to 
groupings  of  data  points  from  the  digitized  profilometer  tracings  (five  data  points  each,  spanning  about 
40  pm  on  the  surface). 
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ratio  HH/HV 
depression  a  <0 


aspect  angle  (degrees) 

Figure  19  G0(HH)/C0(HV)  for  Laboratory  Depression  Angles  of  60° 

Some  of  the  earlier  measurements  obtained  in  the  course  of  this  work  were  made  using  a  rough 
surface  referred  to  as  "original  EBT4"  having  an  average  facet  depth  of  8.68  p.m.  The  slope  angle 
associated  with  this  depth  is  approximately  22°.  At  larger  angles  the  facet  model  deviates  from  the 
measured  values.  Figure  19  is  consistent  with  our  prediction  and  displays  data  for  o0(HH)/o0(HV)  at 
a  60°  laboratory  depression  angle.  This  plot  shows  that  there  is  an  excellent  fit  up  to  a  laboratory 
aspect  angle  of  24°. 

As  a  further  check  on  our  hypothesis  that  curvature  effects  reduce  the  effectiveness  of  the  TP 
approximation,  we  prepared  a  surface  consisting  of  randomly  wound  aluminum  wires  on  an  absorbing 
substrate.  These  measurements  are  displayed  in  Figures  IS  and  16.  Among  the  interesting  features  of 
these  data  is  the  observation  that  at  small  aspect  angles  the  fit  to  the  TP  model  is  good,  and  at  large 
aspect  angles  (Figure  16),  the  data  follow  the  same  pattern  as  the  roughened  surface  at  a  60°  laboratory 
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depression  angle  (Figure  19). 


Figures  11-12  show  plots  of  measurements  for  plastic  in  comparison  with  the  facet  model 
predictions.  The  roughened  plastic  surface  was  prepared  by  pressing  an  acetone  treated  smooth  plastic 
surface  against  a  roughened  aluminum  surface.  After  hardening,  the  plasdc  surface  had  a  roughness 
which  was  verified  by  profilometer  to  be  similar  to  the  rough  metaL  The  complex  index  of  refraction 
of  this  plastic  at  10.6  p.m  was  n  *  1.6  and  k  =  0.002.  The  penetration  depth  is  given  by 

a*1  =  X/4nk  =  422  pm.  This  is  too  large  given  that  the  theoretical  calculations  assume  that  all 
scattering  takes  place  on  the  surface  without  transmission  through  the  facet.  The  scattering  in  the  data 

points  shown  in  Figures  11-12  may  be  due  to  the  transmission  associated  with  small  k. 

To  remedy  this  situation,  we  performed  additional  measurements  on  rough  glass  surfaces. 
Glass  has  a  complex  index  of  refraction  of  n  =  22  and  k  =  0.1  at  10.6  pm.  The  penetration  depth  of 
glass  at  10.6  |im  is  8.4  urn,  sufficiently  small  to  justify  the  neglect  of  transmission  through  facets. 
Two  types  of  rough  glass  surfaces  were  prepared.  Figure  13  displays  scattering  data  from  randomly 
stacked  glass  beads  of  varying  radii  on  a  glass  substrate.  Here  the  deviation  from  the  facet  model  is  at 
very  small  laboratory  aspect  angles.  The  other  set  of  glass  scattering  data  are  displayed  in  Figure  14. 
This  target  surface  was  prepared  by  abrading  a  flat  glass  surface  with  glass  beads  under  pressure. 
Rough  glass  surfaces  prepared  in  this  manner  displayed  many  deep  fissures  when  observed  under  a 
microscope. 

This  surface  did  not  quite  satisfy  the  surface  criteria  necessary  to  test  the  theory.  Nevertheless, 
at  steep  angles  (Figure  14),  the  data  seemed  to  validate  the  facet  model  prediction  up  to  a  laboratory 
aspect  angle  of  15°.  It  is  instructive  to  compare  Figure  14  with  Figure  7  (10/10  data).  Both  appear  to 
have  the  same  range  of  validity.  Both  display  an  HH  null  with  a  depth  of  0.1.  When  comparing 
Figures  13  and  14  with  Figures  7  and  8,  one  observes  that  the  deviation  from  theory  is  more 
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pronounced  for  dielectrics  than  metals  at  larger  laboratory  depression  angles. 

To  obtain  a  better  understanding  of  the  limitations  of  the  TP  approximation,  we  present  some 
numerical  results  of  polarization  state  dependent  scattering  amplitudes  for  both  metal  and  dielectric 
cylinders  of  different  radii. 

The  underlying  rationale  for  doing  this  is  to  juxtapose  the  TP  approximation  with  exact 
calculations  for  curved  surfaces  where  solutions  to  Maxwell's  equations  with  associated  boundary 
conditions  can  be  obtained  numerically.  If  the  tangent  plane  approximation  were  to  be  valid 
independently  of  the  surface  curvature,  one  should  be  able  to  approximate  the  scattering  amplitude 
from  any  cylinder  of  arbitrary  radius  by  the  Fresnel  coefficients  of  the  tangent  plane.  Thus,  in  the  case 
of  a  metal  cylinder,  one  would  conclude  that  the  scattering  intensity  for  a  cylinder  is  independent  of  the 
azimuthal  angle  <j>  and  the  orientation  of  the  incident  linear  polarization.  This,  however,  is  not  the  case. 
An  understanding  of  the  structure  of  the  polarization  dependent  scattering  amplitudes  as  a  function  of 
aA.  is  one  way  to  put  the  tangent  plane  approximation  in  its  proper  perspective. 

The  validity  of  the  tangent  plane  approximation  will  depend  on  how  close  the  ratio  of  HH/VV 
for  a  cylinder  divided  by  HH/W  for  a  plane  is  to  unity.  An  additional  requirement  is  that  the  phase 
differences  for  HH  and  VV  amplitudes  for  planes  and  cylinders  remain  the  same.  It  is  of  interest  to 
note  that  the  HH-W  phase  differences  for  cylinders  and  planes  are  the  same  except  in  the  immediate 
vicinity  of  the  Brewster  angle. 

Figure  20  displays  the  amplitude  ratios  for  metal  cylinders  of  increasing  radii.  The  radius 
ranges  from  3.18  X  to  9.54  X.  The  plots  confirm  qualitatively  the  inequality  of  Eqn.(Al),  ie,  the  larger 
the  radius  of  curvature,  the  better  the  TP  approximation.  One  observes  that  for  large  angles  of 
incidence  the  TP  approximation  fails  even  for  a  radius  of  curvature  on  the  order  of  10  X.  On  the  other 
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hand,  for  steep  angles  of  incidence,  the  TP  approximation  is  effective  even  for  small  radii  of  curvature. 
This  is  consistent  with  our  measurements  on  roughened  aluminum  surfaces,  dielectric  surfaces,  and 
metal  wire  surfaces  at  steep  angles  of  incidence. 

Similar  results  are  displayed  in  Figure  21  for  dielectric  cylinders  with  n  =  1.6  and  k  =  0.1. 
These  differ  from  the  cylinder  plots  for  metals  and  indicate  that  the  TP  approximation  for  dielectrics 
fails  at  smaller  angles  than  for  metals  with  the  same  radius. 

One  notes  that  the  inequality  derived  in  reference  10  doesn’t  distinguish  between  dielectrics  and 

metals. 


Reviewing  our  data  in  light  cf  these  model  calculations,  we  reach  the  following  conclusion. 
The  bulk  of  the  scattered  radiation  ( -  95%  )  is  accounted  for  by  the  facet  model.  There  is,  however,  a 
small  but  finite  component  of  diffuse  radiation  with  a  polarization  that  varies  rapidly  as  a  function  of 
angle.  It  is  the  latter  which  gives  rise  to  deviations  from  the  predictions  of  the  TP  approximation. 
These  deviations  are  especially  pronounced  in  polarization  states  and  at  angles  where  the  TP  model 
predicts  no  scattering  or  small  amounts  of  scattering.  Thus,  at  an  azimuthal  angle  of  0°,  where  the  TP 
model  predicts  a0(HV)  =  0,  a  small  but  finite  <j0(HV)  appears  consistently.  The  other  area  of 
pronounced  deviation  from  the  TP  model  is  at  scattering  angles  that  are  significantly  different  from  the 
specular  direction.  Here,  the  intensity  of  the  scattered  radiation  should  fall  off  exponentially.  This, 
however,  is  not  the  case  due  to  the  underlying  diffuse  background.  The  most  likely  cause  for  this  is 
the  existence  of  areas  on  the  rough  surface  with  small  radii  of  curvature.  This  hypothesis  is  consistent 
with  the  differences  observed  between  metals  and  dielectrics,  and  with  the  scattering  observed  from 
randomly  distributed  wires. 

Other  investigators  have  attributed  depolarization  effects  to  multiple  scattering!7!.  Another 


38 


plausible  source  of  depolarization  may  be  found  in  two  scales  of  roughness!14!.  These  issues  will  be 
addressed  in  a  subsequent  publication. 

We  express  our  thanks  to  M.  Coulombe,  T.  Horgan,  M.  Grund,  M.  Yoon,  E.  Jordan  and 
C  Laramee  for  technical  assistance. 
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Figure  20b  W/HH  for  a  Conducting  Cylinder  with  ka  =  20  divided 
by  W/HH  for  a  Conducting  Plane 
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Figure  20c  VV/HH  for  a  Conducting  Cylinder  with  ka  =  40  divided 
by  W/HH  for  a  Conducting  Plane 
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Figure  20d  VV/HH  for  a  Conducting  Cylinder  with  ka  =  60  divided 
by  W/HH  for  a  Conducting  Plane 
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Figure  21a  W/HH  for  a  Dielectric  Plane 
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Figure  21b  W/HH  for  a  Dielectric  Cylinder  with  ka  *  20  divided 
by  W/HH  for  a  Dielectric  Plane 
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Figure  21c  VV/HH  for  a  Dielectric  Cylinder  with  ka  =  40  divided 
by  VV/HH  for  a  Dielectric  Plane 


Figure  21d  W/HH  for  a  Dielectric  Cylinder  with  ka  =  60  divided 
by  W/HH  for  a  Dielectric  Plane 
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Appendix  A 

An  Algorithm  for  Transforming  Field  Angles  into  Laboratory  Angles 

Having  defined  laboratory  angles  q,  and  V  in  section  II,  we  proceed  with  further  analysis  to 
obtain  an  algorithm  that  will  yield  a  unique  set  of  q,  \|r  for  a  given  set  of  0j,  0s,  an<M>s. 

We  first  show  that  for  a  fixed  0.  and  0$  both  larger  than  0°  (0°  is  a  point  of  degeneracy)  Cos(<j>s) 

can  assume  all  values  between  <1  and  +1.  Denoting  Cos^)  =  x  and 

Cos(0s)  *  y,  equation  (23)  can  be  written  as: 


Cos(<J>s)  = 


-Sin(q)  +  yx 
Vl  -  x2  Vl  -  y2 


(Al) 


For  a  fixed  V  and  £  and  a  fixed  y,  q  is  not  an  independent  variable.  We  therefore  express 
Sin(q)  in  terms  of  x,  y,  and  V/. 


Sin(q)  =  Sin(q-^+^)  =  Sin(q-£)  Cos(£)  +  CosOV^)  Sin(^)  (A2) 

or 


Sin(q)  =  xy  ±  ^Cos2W  "  y2^  '  *2I 

Cos2(y) 


(A3) 


The  minus  sign  in  front  of  the  square  root  has  to  be  inserted  when  (T|-^)  >  90°.  Combining 
(23),  and  (A1)-(A3)  one  writes: 
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Cos(<J>,) 


V7- 


Cos2y 


Vl  -  y2  V 1  -  x2 


(A4) 


We  first  show  that  Cos(<J>s)  can  assume  the  extreme  values  of  ±1.  This  can  be  seen  by  setting 

V  =  0.  The  +  sign  in  front  of  the  square  root  yields  +1,  and  the  -  sign  yields  *1.  It  is  useful  to  rewrite 
the  expression  (14)  as  two  equations  each  valid  in  a  given  range  of  Cos($s).  Assuming  that  y  >  x : 

-xyTanWi'E&'  V77 


Cos((j)f)  > 


Cosfy 
Vl  -y2  Vl  -  x2 


Cosfy 


(A5) 


Eqn.  (A5)  is  valid  for : 


y  Vr?SCos<W£1 


and 


Cos(<})s)  = 


-xyTanV-Y1'^  V^ 


Cosfy 


V 1  -  y2  Vl  -  x2 


(A6) 


Eqn.  (A6)  is  valid  for : 


•1  S  Cos(<J>.)  £  -  i 


s,~  y  V  i-x2 


For  x  >  y,  -  ^  ~~^2  must  be  replaced  with  - ^ .  Since  both  functions  are 
continuous  functions  of  y,  for  0  £  y  £  Cos'^y),  Cos(<J>s)  assumes  all  values  in  between. 

The  algorithm  for  determining  q,  y  for  a  given  9.,^,  ^  is  as  foUows.  Compute  Cos(4>s) 
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and  select  the  larger  of  Cos(6s)  or  Cos(8j),  thereby  setting  the  appropriate  range  for  equations  (AS) 
and  (A6).  The  value  of  Cos(<J>s)  determines  whether  Eqn.  (A5)  or  (A6)  is  to  be  used  Having  selected 
the  correct  equation,  one  solves  for  Cos(y)  and  determines  \jr.  Given  Cos(y),  one  then  can  determine 
Cos(^)  or  rj  remains  to  be  determined  and  is  obtained  from:  y  =  Cos (0,)  =  SinCn-£)  Cos(y). 

Since  Sin(r| — is  symmetric  around  90° ,  a  unique  determination  of  t|  depends  on  whether 
Eqn.  (A5)  or  (A 6)  was  used  For  Eqn.  (A5),  £  90°  and  for  Eqn.  (A6)  <,  90°. 

To  complete  this  section,  we  write  down  the  direction  of  linear  polarization  that  corresponds  to 
an  H  transmit  state.  It  is  given  by: 

r,  x  n"  =  -ii  Sin(v)  +  »2  Sin(^)  Cos(y)  (A7) 

An  H  receive  state  is  given  by: 

I'D  *  n"  =  -[ii  Sin(-n)  Cos(^)  +  b  Cos(tj-^)]  Cos(y)  + 13  Cos(ri)  Sin(y)  (A8) 

The  V  receive  state  is  obtained  from: 

rDx(rDxn^  =  rD(rD-  «0-n"  (A9) 


1*0  x  (rD  x  n")  =  ii  [Cos(n)  Sindi-^)  +  Sin(5)]  Cos(y)  +  i2£Sin(T|)  Sin(T|-^>  *-  (A10) 

Cos(v)  -  Sin(\|/)]  -  i3  Cos(^)  Cos(y) 
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Appendix  B 

Additional  Scattering  Data 

To  conoborate  our  hypothesis  that  the  deviation  from  the  TP  A  (Tangent  Plane  Approximation) 
model  is  due  to  the  small  amounts  of  surface  curvature,  we  prepared  a  set  of  roughened  aluminum 
surfaces  that  were  blasted  with  different  size  glass  beads.  Half  of  the  surfaces  were  electropolished 
after  blasting,  the  other  half  were  left  untouched  after  blasting. 

Figures  AC1  (Appendix  Cl)  to  Figure  AC14  display  data  in  laboratory  angles  as  a  function  of 
laboratory  aspect  angle.  The  roughened  surfaces  from  which  these  data  were  obtained  is  designated 
EBT4.  E  stands  for  electropolished  and  BT4  signifies  that  the  surface  was  blasted  with  glass  beads 
BT4. 

Figures  AC15  to  AC26  display  data  from  the  same  surface  in  field  angles.  Here  the  incident 
and  scattered  polar  angles  are  fixed  and  the  azimuthal  scattering  angle  is  varied. 

Figures  AC27  to  AC40  contain  plots  as  a  function  of  laboratory  aspect  angles  from  another 
surface  designated  by  EBT6.  Figures  AC41  to  AC54  represent  data  from  EBT6  in  field  angles. 
Figures  AC55  to  AC68  display  data  in  laboratory  coordinates  from  a  third  surface,  designated  EBT13. 
Figures  AC69  to  AC82  contain  scattering  data  in  field  coordinates  for  surface  EBT  4+13. 

In  figures  AC109  to  AC162  scattering  data  for  laboratory  angles  from  the  four  different 
surfaces  NBT4,  NBT6,  NBT13,  and  NBT4+13  are  displayed.  N  stands  for  non-electropolished.  The 
designation  BT4,  BT6,  BT13,  and  BT4+13  refer  to  three  different  glass  beads  used  in  the  sand 
blasting  process. 

Figures  AC163  to  170  are  histograms  displaying  percent  of  surface  Vs  slope. 

Tables  I  and  II  summarize  profilometric  surface  data  including  the  percentage  of  surface  that 
has  a  curvature  of  less  than  3X  (31.8|im).  Tables  HI  and  IV  display  the  depth  of  the  HH  nulls.  The 
TPA  predicts  an  absolute  zero  at  the  location  of  the  null,  and  in  practice  would  be  limited  by 
instrumental  uncertainty,  ie.,  the  extinction  ratio  of  the  polarization  state  analyzer.  If  the  measured  null 
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is  above  the  instrumental  threshold,  it  would  indicate  that  another  source  is  present  In  Tables  II  and 
IQ  we  demonstrate  another  correlation  between  the  depth  of  die  HH  nulls  and  percentage  of  curvature. 
The  smaller  the  percentage  of  curvature,  ( r  <  3X ),  the  deeper  the  HH  null. 
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AZIMUTHAL  ANGLE  (DEGREES) 

FIG.  AC-25 
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Intensity 


a0(HH)/o0(HV)  for  laboratory  depression  angles  =  2;s  =  30 


a0(HH)  and  c0(HV)  for  laboratory  depression  angles  ^.  =  ^  =  30 
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a0(HH)/a0(HV)  for  laboratory  depression  angles  ^  =  40 


<70(HH)  and  c0(HV)  for  laboratory  depression  angles  =  £  =  40 
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V**V*.. 


HH/HV  data 
HH/HV  Model 


AC-38 


Aspect  angle  (degrees) 


a0(HH)  and  a0(HV)  for  laboratory  depression  angles  £,.  =  =  50 


lC-39 


Aspect  angle  (degrees) 


a0(HH)/a0(HV)  for  laboratory  depression  angles  =  60 


o0(HH)  and  c0(HV)  for  laboratory  depression  angles  ^  =60 


10  DEGREE  POLAR  ANGLE 
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a0(HH)/o0(HV)  for  laboratory  depression  angles  =  50 
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a0(HH)  and  cr0(HV)  for  laboratory  depression  angles  =  50 
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o0(HH)/ao(HV)  for  laboratory  depression  angles  =  £  =60 
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o0(HH)/a0(HV)  for  laboratory  depression  angles  !;.  =  £,  =20 
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o0(HH)/oo(HV)  for  laboratory  depression  angles  ^  =  4S  =  30 


30/30  HH  and  HV  for  NBT  4 
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g0(HH)  and  g0(HV)  for  laboratory  depression  angles  =  £  =  10 
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TABLE  I 


*  Radius  oi  Curvalura 
%  of  Surfaca  <  6  lambda 
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limensional  tracings  in  three  directions  through  a  point  on  the  surface;  %  means  percent  of  total  number  of  data  points;  lambda  =  10.6  microns. 


Table  II 

Distribution  of  Slopes  for  Electropolished  and  Non-Electropoiished  Surfaces 


*  Surface  slope 
%  of  Surface  <  0.25 

95.6 

97.5 

98.9 

95.8 

99.8 

99.6 

97.6 

96.8 

*  Surface  slope 
%  of  Surface  <  0.20 

91.8 

94.9 

97.1 

93.6 

99.7 

98.9 

95.4 

93.7 

*  Surface  slope 
%  of  Surface  <0.15 

82.4 

88.5 

93.2 

87.9 

98.2 

96.7 

90.7 

85.8 

*  Surface  slope 
%  of  Surface  <  0..10 

68.8 

73.5 

ZZL 

*  28 

89.5 

88.6 

79.5 

72.2 

*  Surface  slope 
%  of  Surface  <  0.05 

40.1 

45.5 

49.2 

46.0 

9  68 

9  99 

52.8 

42.9 

STANDARD 

DEVIATION 

0.126 

0.09 

860  0 

690  0 

0.047 

0.055 

980  0 

980  0 

RMS 

SLOPE 

0.169 

0.129 

0.102 

0.134 

0  07 

0.08 

0.118 

0  125 

■ 

0.112 

0.091 

0.075 

0.094 

0.052 

0.058 

0.081 

0.092 

GLASS  BEAD  (DIA.) 

E  BT4  (500) 

N  BT4  (500) 

E  BT6  (250) 

N  BT6  (250) 

E  BT13  (50) 

N  BT13  (50) 

X  X 

1  5 

«  n 

H  j" 

m  is 

♦  + 

O  O 

U1  Z 
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Table  IH  Null  Depth  at  Three  angles  for  Surfaces  Studied 


Incident/Scattered  anales 

EBT4 

EBT6 

EBT13 

EBT4+13 

EBT4  (original) 

7.517.5 

0.015 

0.01 

0.007 

0.009 

0.026 

10/10. 

0.015 

0.006 

0.007 

0.011 

0.021 

20/20 

0.054 

0.02 

0.008 

0.035 

0.027 

SUMOFAVG. 

0.084 

0.036 

0.022 

0.055 

0.074 

Table  IV  Surfaces  Ranked  from  Deepest  to  Shallowest 

by  Depth  of  Null 
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OF 

ROME  LABORATORY 

Rome  Laboratory  plans  and  executes  an  interdisciplinary  program  in  re¬ 
search,  development ,  test ,  and  technology  transition  in  support  of  Air 
Force  Command,  Control,  Communications  and  Intelligence  (C  I)  activities 
for  all  Air  Force  platforms.  It  also  executes  selected  acquisition  programs 
in  several  areas  of  expertise.  Technical  and  engineering  support  within 
areas  of  competence  is  provided  to  ESD  Program  Offices  (POs)  and  other 
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ESD  elements  to  perform  effective  acquisition  of  C  l  systems.  In  addition, 
Rome  Laboratory's  technology  supports  other  AFSC  Product  Divisions,  the 
Air  Force  user  community,  and  other  DOD  and  non-DO D  agencies.  Rome 
Laboratory  maintains  technical  competence  and  research  programs  in  areas 
including,  but  not  limited  to,  communications,  command  and  control,  battle 
management,  intelligence  information  processing,  computational  sciences 
and  software  producibility,  wide  area  surveillance/sensors,  signal  proces¬ 
sing,  solid  state  sciences,  photonics,  electromagnetic  technology,  super¬ 
conductivity,  and  electronic  reliability/maintainabUity  and  testability. 


